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1. Introduction 

To give a possible explanation of color confinement G. 't Hooft [1] and S. Mandelstam [2] 
suggested long time ago that vacuum in gauge theories behaves like a magnetic (dual) super- 
conductor. The dual superconductivity hypothesis relies upon the very general assumption 
that the dual superconductivity of the ground state is realized if there is condensation of 
Abelian magnetic monopoles. Indeed, lattice calculations have given evidence of Abelian 
magnetic monopoles condensation [3-13]. 

On the other hand numerical evidence has emerged [14-27] in favor of the so called 
center vortex picture where the vacuum consists of a coherent condensate of magnetic flux 
tubes. Also this theoretical proposal has been advanced long time ago [28-35]. 

In this paper we compare the dual superconductivity scenario with the vortex con- 
densate picture both in Abelian and non Abelian pure gauge lattice theories. A partial 
account of results discussed in the present paper has been reported in Ref. [36]. In this 
work we do not consider the center vortices. Instead we study the Abelian vortices which 
eventually give rise to a coherent condensate of Abelian magnetic flux tubes. 

Monopole or vortex condensation can be detected using order /disorder parameters [37- 
39]. At zero temperature we employ a disorder parameter defined in terms of a gauge- 
invariant lattice effective action introduced [40-42] using the lattice Schrodinger func- 
tional [43-46]. At finite temperature the disorder parameter has been defined [7] in terms 
of a thermal partition functional. 

For the sake of clarity let us recall our implementation of the gauge-invariant lattice 
Schrodinger functional previously discussed in Ref. [42]. Let us consider the continuum 
Euclidean Schrodinger functional in Yang-Mills theories without matter field: 

Z[A(/),A«] = (a(/) |e-^^P| A«) . (1.1) 
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In Eq. (|r3) H is the pure gauge Yang-Mills Hamiltonian in the temporal gauge, T is the 
Euclidean time extension, while V projects onto the physical states. and A'^^'^\x) 

are static classical gauge fields, and the state |A) is such that 

(A|'f) = *[A]. (1.2) 

Inserting an orthonormal basis \^n) of gauge invariant energy eigenstates in Eq. ( |1.1D 

^[A(^), A«] = ^e-^"^^'„[A(^)]**[A»] . (1.3) 

n 

Since we are interested in the lattice version of the Schrodinger functional, it makes sense 
to perform a discrete sum in Eq. ( |1.3| ) for the spectrum is discrete in a finite volume. 
Eq. (1.3) shows that the Schrodinger functional is invariant under arbitrary static gauge 
transformations of the fields A^-^^ and A^*^ 

Using standard formal manipulations and the gauge invariance of the Schrodinger 
functional it is easy to rewrite Z[A^-f \ A^'^^] as a functional integral [47-49] 

with the constraints: 

A(xo = 0) = aW , 



A(2;o = T) = A(^) 



(1.5) 



Strictly speaking we should include in Eq. ( |1.4D the sum over topological inequivalent 
classes. However, it turns out that [43-46] on the lattice such an average is not needed 
because the functional integral Eq. (1^) is already invariant under arbitrary gauge trans- 
formations of A^*) and A^-^^ 

On the lattice the natural relation between the continuum gauge fields and the corre- 
sponding lattice links is given by 

Ufj_{x) = F ex.p i^iag J dt Af^{x + atfl)^ (1-6) 

where P is the path-ordering operator, a is the lattice spacing and g the bare gauge coupling 
constant. 

The lattice implementation [43-46] of the Schrodinger functional, Eq. ( |1.4| ), is now 
straightforward: 

Z[U^-f\U^ = J VU^ e-^. (1.7) 

In Eq. ( |l.7D the functional integration is done over the links U^{x) with the fixed boundary 
values: 

Uk{x)\x,=o = u'^\ Uk{x)l,=T = u\f^ ■ (fc = l,2,3). (1.8) 
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Links in temporal direction are not constrained. S is the standard Wilson action modi- 
fied [43-46] to take into account the boundaries at X4 = 0, T. For SU(N) 

S = PY. ^/^-(^) {1 - ^I^eTr[l - U^u{x)]\ , /? = ^ , (1.9) 

where Ufj,u{x) are the plaquettes in the (/i, i/)-plane and the weights Wf^uix) are given by 

, , I 1/2 spatial plaquettes at X4 = 0, T ,^ ^. 

Wf,u[x) = ( _ . (1.10) 

I 1 otherwise 

It is possible to improve the lattice action S by modifying the weights W^j/'s [43-46]. Note 
that, due to the fact that / U^f\ one cannot impose periodic boundary conditions 
in the Euclidean time direction. On the other hand one can assume periodic boundary 
conditions in the spatial directions. 

Let us consider, now, a static external background field A°^*(x) = A^^*(x)Aa/2, where 
Aa/2 are the generators of the SU(N) Lie algebra. We introduce a new functional: 

r[A'=^*] = -iln|^|^\ , (1.11) 



T { Z[0] 
where 

Z[C/^^*] = Z[[/^^*, [/^^*] , (1.12) 

and Z[0] is the Schrodinger functional Eq. (|1.12|) with A°^* = (C/^^* = 1). The lattice 
link C7|^* is obtained from the continuum background field A°^* through Eq. (|1.6| ). 

From the previous discussion it is clear that r[A''^*] is invariant for lattice gauge 
transformations of the external links U^^^. Moreover, from Eq. ( |1.3| ) it follows that 

lim r[A^^*] = ^o[A"=^*] - ^o[0] (1.13) 

T — >oo 

where £'o[A'^^*] is the vacuum energy in presence of the external background field. In other 
words r[A'^^*] is the lattice gauge-invariant effective action for the static background field 

^cxt 

Note that, since our definition of the lattice effective action uses the lattice Schrodinger 
functional with the same boundary fields at X4 = and X4 = T, we can glue the two 
hyperplanes X4 = and X4 = T together. This way we end up in a lattice with periodic 
boundary conditions in time direction too. Therefore our lattice Schrodinger functional 
turns out to be 

Z[Un = I VU, e-^ , (1.14) 

where the functional integral is defined over a four-dimensional hypertorus with the "cold- 
wall" 

Uk{xi,x%,=^ = Ul^\x). (1.15) 
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Moreover, due to the lacking of free boundaries, the lattice action in Eq. ( |1.14| ) is now the 
familiar Wilson action 



(1.16) 



We also impose that links at the spatial boundaries (we mean the spatial boundaries of 
each given time slice and not only spatial boundaries of the slice X4 = 0) are fixed according 
to Eq. (1.15). In the continuum this last condition amounts to the usual requirement that 
the fluctuations over the background field vanish at infinity. 



We want to extend our definition of lattice effective action to gauge systems at finite 
temperature. In this case the relevant quantity is the thermal partition function. In the 
continuum we have: 



Tr 



PA (A 



(1.17) 



where (3t is the inverse of the physical temperature, H is the Hamiltonian, and V projects 
onto the physical states. As is well known, the thermal partition function can be written 
as [49]: 



Tr 



VA^{xi, x) e- ^^^4 / d^xCY-Mix,xi) 



Afj,{l3T,x)=Afj,{0,x) 



On the lattice we have: 



Tr 



-I3tH 



IU^{I3t,x)=U^{0,x) 
Comparing Eq. (|1.19|) with Eqs. ( 1.14 ) and (|1.15|) , we get: 

-I3tH _ 



PC/^(x4,^)e-^^. 



(1.18) 



(1.19) 



Tr 



VUk{x) Z[Uk{x)] , 



(1.20) 



where Z[Uk{x)] is the Schrodinger functional Eq. ( |1.14 ) defined on a lattice with L4 = /3t, 
with "external" links Uk{x) at X4 = 0. We recall that, by definition, Z[Uk{x)] includes the 
integration over U4{0,x). 

We are interested in the thermal partition function in presence of a given static back- 
ground field A^^^{x). In the continuum this can be obtained by splitting the gauge field 
into the background field A^^*(x) and the fiuctuating fields r]{x). So that we could write 
formally for the thermal partition function 2'jn[A'^^*]: 



(1.21) 



Actually, to give a meaning to Eq. (|1.21 ) one must define the states | A*^^* ,77 > . For 
instance, in a perturbative approach, a meaningful definition of 2^t[A'^^*] is given by the 
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background field method. In the non-perturbative lattice approach the implementation of 
Eq. (1.21) could be obtained as follows. We write 



(1.22) 



where U^^{x) is related to A'^^*(x) by Eq. (|T]B) and the f/fc(x)'s are the fluctuating links 



(related to the fluctuating fields ri{x)) . Thus Eq. ( 1.20[ ) suggests the following definition 

Zr[A-*] = / VUk{x) Z[^r*(x),t/fc(x)] , (1.23) 



where we integrate over the fluctuating links Uk{x), while the U^'^{x) links are flxed. Note 



that in Eq. ( 1.23| ) only the spatial links exiting from sites belonging to the hyperplane X4 = 
are written as products of the external links U^^{x) and the fluctuating links Uk{x). They 
will be called "frozen links", while the remainder will be called "dynamical links". From 
the physical point of view we are considering the gauge system at finite temperature in 
interaction with a fixed external background field. Therefore in the Wilson action Sw we 
do not include the contribution of plaquettes built up with only frozen links. The temporal 
links C/4(x4 = 0, x) are not constrained and satisfy the usual periodic boundary conditions. 
We stress that p.b.c.'s in temporal direction are crucial to retain the physical interpretation 
of the functional ^^^'[A'^^*'] as thermal partition function. 
Now, it is easy to see that in Eq. ( |1.23| ) we have 



ui-\x)Mx) =z[ui-\. 



(1.24) 



Indeed, let us consider an arbitrary frozen link U^'^{x)Uk{x). This link enters in the Wilson 
action through the plaquette: 

PUxi = 0,f) = Tr{;7r*(x)^fc(x)C/4(0,x + ^)C/2(l,x)C/i(0,x)} . (1.25) 



Now we observe that the link [74(0, x + k) in Eq. ( 1.25| ) is a dynamical one, i.e. we are 
integrating over it. So that the dependence on Uk{x) can be re-absorbed by a change of 
integration variable. Therefore we obtain 

Pkiixi = 0,x) = Ti\ui?'\x)Ui{Q,x + k)ul{l,x)ul{Q,x)^ . (1.26) 



It is evident that Eq. (1.26) in turns implies Eq. ( [L.24 ). Then, we see that in Eq. (jO^ ) 
the integration over the fiuctuating links U{x) gives an irrelevant multiplicative constant. 
So that we are led to define the lattice thermal partition function in presence of a static 
background field as 



Zt [a-*] = [ 
Jv 



VUe 



Uk{PT.S)=Uk{0,x)=Ul^\x) 



(1.27) 



Note that the thermal partition function 2^t[A'^^*], as defined by Eq. ( 1.27] ), is invariant 
for time-independent gauge transformations of the background field A*^^*. On a lattice 
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with finite spatial extension we impose that the spatial links exiting from sites belonging 
to the spatial boundary of a generic time slice X4 {x4 7^ 0) are fixed according to the 
boundary conditions Eq. (|1.15| ) used in the Schrodinger functional. Thus we see that, 
sending physical temperature to zero, the thermal functional Eq. ( 1.27| ) reduces to the 
zero-temperature Schrodinger functional Eq. ( [L.14 ) on a finite lattice. 



Let us now introduce our disorder parameter to detect monopole or vortex condensa- 
tion. At zero-temperature [7] 

^ = e-^i' f-^^ = , (1.28) 

where 2^[A''^*] is the lattice Schrodinger functional with monopole or vortex static back- 
ground field A*^^* . According to the physical interpretation of the effective action Eq. ( |1.11| ) 
-Eb.f. is the energy to create a monopole or a vortex in the quantum vacuum. If there is 
condensation, then E'b.f. = and fj, = 1. 

At finite temperature the disorder parameter turns out to be related to the monopole or 
vortex free energy [7,36]. In particular, the finite temperature disorder parameter is defined 
by means of the thermal partition function in presence of the given static background field 



Eq. (11:271) : 



^ = e-Fb.f./W = '^^^^^ ^ , (1.29) 



From Eq. (|]29|) it is now clear that Fy^j^ is the free energy to create a monopole or a vortex. 



If there is condensation, then i<b,f, = and fi = 1. 

As already stated, our disorder parameter is gauge-invariant for time-independent 
gauge transformations of the external background fields, since it has been defined in terms 
of the Schrodinger functional. Let us stress that gauge invariance for time-independent 
gauge transformations of the external background fields implies that we do not need to do 
any gauge fixing to perform the Abelian projection. Indeed, after choosing the Abelian 
direction, needed to define the Abelian monopole or vortex fields through the Abelian 
projection, due to gauge invariance of Schrodinger functional for transformations of back- 
ground field, our results do not depend on the selected Abelian direction, which, actually, 
can be varied by a gauge transformation. 

The plan of the paper is the following. In Sect. 2 we study the condensation of vortices 
and monopoles in the zero temperature lattice U(l) pure gauge theory. In Sect. 3 we 
compare Abelian monopole and vortex condensation for finite temperature SU(2) lattice 
gauge theory. Sect. 4 is devoted to finite temperature SU(3) gauge theory where, according 
to the choice of the Abelian subgroup, two different kinds of Abelian monopoles and vortices 
can be defined. Our conclusions are drawn in Sect. 5. 

2. U(l) 

In this Section we study the monopole and vortex condensation in lattice pure gauge U(l) 
theory at zero physical temperature. The disorder parameter Eq. ( |1.28 ) is defined in terms 
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Figure 1: The derivative of the energy to create a monopole (nmon = 1), Eq. versus /3 for 

U(l) lattice gauge theory on a L"^ lattice (circles refer to L = 16, squares to L = 24, and diamonds 
to i 32). 



of the lattice effective action Eq. ( 1.11| ) with a Dirac magnetic monopole or a magnetic 
vortex background field. 

Let us start with the case of the Dirac magnetic monopole background field. In the 
continuum the magnetic monopole field with the Dirac string in the direction n is 



eb{f) 



r X n 



2 r(r 



(2.1) 



r ■ n] 



where, according to the Dirac quantization condition, nmon is an integer and e is the electric 
charge (magnetic charge = nmon/2e). We consider the gauge-invariant background field 
action Eq. ( |1.11| ) where the external background field is given by the lattice version of the 
Dirac magnetic monopole field. By choosing n = X3 we get: 



C/f|(f) = cos[05^^°(f)] + isin[e5^°"(^)] 



(2.2) 
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with 

{X2 - X2) 1 



2 l^^monl l^monl (^^3 ^2 

^mon^^-j _ I ^mon (^1 ~ ^l) 1 

2 l^^monl I -^1x1011 1 (•^S 

In Equation ( |2.3| ) (^i, A^2)^3) s-re the monopole coordinates and Xmon = (a; — X). In the 
numerical simulations we put the lattice Dirac monopole at the center of the time slice 
X4 = 0. To avoid the singularity due to the Dirac string we locate the monopole between 
two neighboring sites. We have checked that the numerical results are not too sensitive to 
the precise position of the magnetic monopole. 

To avoid the problem of dealing with a partition function we consider -Emon = dEmon/df]. 
It is easy to see that E'^^^ is given by the difference between the average plaquette < U^i, > 
obtained in turn from configurations without and with the monopole field: 

-E'mon = y[< Ufj,u >nmon=0 " < Uf_,u >nmonf^o] > (2-4) 

where V is the spatial volume. 

We performed lattice simulations on 16^, 24^ and 32^ lattices using an APEIOO com- 
puter. The spatial links belonging to the time slice X4 = and to the spatial boundaries 
of the other time slices (X4 7^ 0) are constrained, therefore they are not updated dur- 
ing Monte Carlo simulation. Therefore we must consider only the "dynamical links" in 
the computation of E'^^^. This means that the generic plaquette U^y{x) = Ufj_{x)Uy{x + 
fl)Uji{x + i')Ul{x) contributes to Eq. (|2.4|) if the link U^{x) is a "dynamical" one (i.e. it is 
not constrained in the lattice Schrodinger functional integration). 

Since we are measuring a local quantity such as average plaquette < U^y > , a low statistics 
(from 1000 up to 5000 configurations) is required in order to get a good estimate of E'^^^. 
Statistical errors have been estimated using the jackknife resampling method [50,51], mod- 
ified to take into account the statistical correlations between lattice configurations. 
In Fig. 1 we report our numerical results for E'^^^ versus P for three different lattice sizes. 
We see that in strong coupling region (/? < 1), monopole internal energy derivative is 



zero, insensitive to the lattice size. This means that, according to Eq. (1.28), the disorder 
parameter /i ~ 1. On the other hand, near the critical coupling /3c ^ 1, £'mon displays a 
sharp peak which increases by increasing the lattice volume. In the weak coupling region 
(/3 Pc) the plateau in E'^^^ indicates that the monopole energy tends to the classical 
monopole action which behaves linearly in /3. 

In order to obtain fi we perform the numerical integration of E'^^^^ 

Emon = / ^mon '^P (2-5) 







In Fig. 2 we display the logarithm of the disorder parameter ^ versus /3. In the confined 



phase we see that In/i = 0, so that Eq. ( 1.28 ) tells us that the energy required to create a 



monopole is zero and therefore monopoles condense in the confined U(l) vacuum. Corre- 
spondingly the disorder parameter ^ is different from zero in the confined phase. Moreover 
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Figure 2: The logarithm of the disorder parameter, Eq. ( 1.28 ), versus /3 for U(l) lattice gauge 
theory on a L'^ lattice. Symbols as in Fig. 1. 



numerical data suggest that /i — > when /3 ^ /3c in the thermodynamic limit. Note that 
the different curves for fi corresponding to increasing lattice sizes seem to cross, suggesting 
a first order phase transition. However, in order to extract the critical parameters and to 
determine the order of transition we need to perform a finite size scaling analysis, which 
will be addressed in a future work. 

Let us now consider the case of U(l) magnetic vortices. The continuum gauge potential 
for a classical magnetic vortex along the xs-direction with nvort units of elementary flux 
(/) = 27r/e is given by 



' ~ e (Xi)2 + (X2)2' 
^ext ^ xi .2.6) 

2 e (Xi)2 + (X2)2' 

Af^ = . 
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Figure 3: The derivative of the energy to create a monopole or a vortex versus fi for U(l) lattice 
gauge theory on a 16^ lattice. Open circles refer to a monopole background field (rimon = 1) and 
full circles to a vortex background field (nvort = 1) • 



The corresponding lattice links are: 





= ^O^K2 


/]Vort 

£^1,2 


— "F '^vort ^ 


ur\x) 


= 1. 



X2A 



(2.7) 



As in the monopole case we evaluated numerically the /3-derivative, -Eyorf energy 
to create a vortex. In Fig. 3 we compare the monopole and vortex energy derivative for 
a 16^ lattice . While we clearly see that in the strong coupling monopoles condense, the 
/3-derivative of the vortex energy displays rapid oscillations whose amplitude seems to 
increase near the critical coupling. This means that we cannot obtain a reliable estimate 
of ii'vort by a numerical integration of E'^^^^. However, we can safely say that our data do 
not show any signal of vortex condensation which would imply that E'^^n = at strong 
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Figure 4: The derivative of the free energy, Eq. (3.5), versus P for monopoles (open circles) and 
vortices (full circles) for SU(2) on a 24^ x 4 lattice. The absolute value of the Polyakov loop. 



Eq. (3.7), is also displayed (open squares). 



couplings. 

Thus, we may conclude that in U(l) lattice theory the strong coupling confined phase is 
intimately related to magnetic monopole condensation. 



3. SU(2) 

In a previous work [7] we studied Abelian magnetic monopole condensation in finite tem- 
perature SU(2) lattice gauge theory. For SU(2) the maximal Abelian group is an Abelian 
U(l) group. In the continuum the Abelian monopole field is given by 



2 \x\ (\x\ — X ■ n) 



(3.1) 



where ft is the direction of the Dirac string and, according to the Dirac quantization 
condition, integer. The lattice links corresponding to the Abelian monopole 
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field Eq. ( |3.lD can be readily obtained from Eq. (|l.q ). By choosing n = X3 we have: 

u!f{x) = cos[eT,n^)] + icT3sin[e?;r(^)] , 

ur\x) = i, 



with 



imon, 



(3.2) 

{X2 - X2) 1 
4 |2^mon| l^monl (^^3 ^s) 

^^■^ _ I "-mon (a;i - Xi) 1 

4 l^^monl l^monl (^^3 ^s) 

where (Xi,X2,X3) are the monopole coordinates, Xmon = (x — X) and the Ua's are the 
Pauli matrices. 

As discussed in Sect. 1, at finite temperature the disorder parameter, Eq. ( 1.2g| ), is 
defined by means of the thermal partition function in presence of the Abelian 

monopole background field Eq. (3.1). Numerical results in Ref. [7] show that the monopole 
disorder parameter ^ is different from zero in the confined phase and suggest that it tends 
to zero when approaching the critical coupling in the thermodynamic limit. Thus SU(2) 
confining vacuum does display the Abelian monopole condensation in accordance with dual 
superconductivity hypothesis. 

Now we want discuss what happens if we consider an Abelian vortex background field. 
In SU(2) gauge theory, Abelian vortex field on the lattice is given by 

Uffix) = cosl^r^'H^')] + ^^3 MOlfix)] , 



"^'^ '-^'-^ 2 (Xi)2 + (X2)2 ■ 

Again, the derivative of the free energy required to create a vortex 

/ d 

Kovt = -Q^Fmrt , (3.5) 

can be easily evaluated as the difference between the average plaquette < U^y > without 
the vortex background field (i.e. nvort = 0) and with the vortex background field (nvort 7^ 0) 

^vort = V[< U^,u >nvon=0 " < Uf,^ >n_t^o] , (3.6) 

where V is the spatial volume. In Eq. ( |3.6| ) we include only the contributions due to the 
dynamical links (see discussion after Eq. (p.4|)). A low statistics (from 2000 up to 10000 
configurations) is required in order to get a good estimate of the derivative of the free 
energy. In Fig. 4 we display the derivative of the vortex free energy versus /3 for nvort = 10 
on lattices with = 4 and Ls = 24. We see that -F^ort vanishes at strong coupling and 
displays a rather sharp peak near (3 2.2. We expect that this peak corresponds to the 
finite temperature deconfinement transition. In Fig. 4 we also display the absolute value 
of the Polyakov loop in time direction 

^ = ^E^Tr n U,{x,,x) (3.7) 

X X4 = l 
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Figure 5: The logarithm of the disorder parameter, Eq. ( 1.29 ), versus /3 for vortices (full circles) 
and monopoles (open circles). Data refer to SU(2) gauge theory on a x 4 lattice. 



and, indeed, we can see that the peak in F^^^.^ corresponds to the rise of the Polyakov loop. 
In weak couphng region the plateau in i^voj-t indicates that vortex free energy tends to the 
classical vortex action which behaves linearly in f3. In Fig. 4 we display for comparison 
also the derivative of the Abelian monopole free energy for nmon = 10. It turns out that 
monopoles and vortices data agree within statistical errors. To appreciate better this last 
point we plotted in Fig. 5 the logarithm of the disorder parameter Eq. ( 1.29| ) for monopoles 
and vortices respectively. Fig. 5 shows clearly that monopoles and vortices agree quite 
perfectly. Our results strongly suggest that also Abelian vortices could play a role in the 
dynamics of confinement. 



4. SU(3) 



For SU(3) gauge theory the maximal Abelian group is U(l)xU(l), therefore we may intro- 
duce two independent types of Abelian monopoles or Abelian vortices. 
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Figure 6: The derivative of the free energy versus f3 for monopoles (open circles) and vortices (full 



circles). The absolute value of the Polyakov loop, Eq. (15), is also displayed (open squares). Data 
refer to SU(3) gauge theory on a 32'^ x 4 lattice. 



Let us consider the Abelian monopole field. The first type of Abelian monopole field is 
derived considering the A3 diagonal generator, we name it Ts Abelian monopole (following 
Ref. [7]). On the lattice it is given by 



e-*^?:r(-) 

1 



(4.1) 



with 6*5^1'^ (f) defined in Eq. (pl) . 

The second type of independent Abelian monopole can be obtained by considering the 
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Figure 7: The logarithm of the disorder parameter, Eq. ( 1.29 ), versus (3 for Ts monopoles (open 
circles) and Tg vortices (full circles). Data refer to SU(3) gauge theory on a 32'^ x 4 lattice. 



diagonal generator Ag. In this case we have the Tg Abelian monopole: 



Utix) = 1 



e 




e 







with 
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71 
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4 l^^monl l^^mon I - (X3 - ^3) 

(xi - Xi) 1 



moni l-^mon 



^monl - {X3 " ^s) 



(4.2) 



(4.3) 
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Figure 8: The free energy derivative for T3 vortices (open circles) and Tg vortices (full circles) 
versus /?. Data refer to SU(3) gauge theory on a 32^ x 4 lattice. 



Analogously, we have the Ts Abelian vortex: 



e-*^™"(^^) 
1 



J7|^*(f) = 1 



(4.4) 



2 (Xi)2 + (X2)2 
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Figure 9: The free energy derivative for Tg vortices (open circles) and Tg vortices (full circles) 
versus /?. Data refer to SU(3) gauge theory on a 32^ x 4 lattice. 

and the Tg Abelian vortex: 















e' 


















QVort 



U^ix) = 1 



(4.5) 



IT = T 



1 ^vort X2,l 

V3 2 (xi)2 + (a;2)2 ' 



Other Abelian monopoles and vortices can be generated by considering linear combinations 
of A3 and As generators. In particular we have also considered the following two linear 
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Figure 10: The free energy derivative for T3 vortices (open circles) and T^a vortices (full circles) 
versus f3. Data refer to SU(3) gauge theory on a 32'^ x 4 lattice. 



combinations of A3/2 and A8/2 



2 2 2 2 



and 



^ 2 2 2 2 
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oi 
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1 

2V3- 



(4.6) 



(4.7) 



The linear combinations given in Eq. ( [4 .61 ) and in Eq. (4/?) have been also considered in 
Ref. [4,5] and in Ref. [25] respectively. 



-18- 



200 



150 



c 
o 
V E 



100 



50 







I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



O Tg 



SU(3) 



123456789 10 

P 

Figure 11: The free energy derivative for Tg monopoles (open circles) and monopoles (full 
circles) versus /3. Data refer to SU(3) gauge theory on a 32^ x 4 lattice. 



For reader convenience, let us summarize the main results of our previous investigation 
of Abelian monopole condensation [7] in SU(3). In Ref. [7] we found that there is condensa- 
tion of AbeUan monopoles in the non-perturbative vacuum and that SU(3) vacuum reacts 
moderately strongly in the case of the Tg Abelian monopole (in particular the peak value 
of -Fmon foj" ^8 Abelian monopoles is about two times higher than the corresponding peak 
for Ts Abelian monopoles). 

In the present paper we compare Tg Abelian monopoles and vortices. In Fig. 6 the 
free energy derivative for monopoles (with rimon = 10) and vortices (with nvort = 10) is 
displayed versus /3 for a lattice with Lg = 32 and Lt = 4. We also display the absolute 
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value of the Polyakov loop in the time direction 



^ = ^E^'^ n U,{x,,x). (4.8) 

X X4 = l 

As can be argued from Fig. 6, F^^^.^ behaves like Fl^^^. Indeed, the free energy derivatives 
are zero within errors in the strong coupling region and display a sharp peak in corre- 
spondence of the rise of the Polyakov loop. In the weak coupling region the free energy 
derivatives are almost constant. The values of the plateau correspond to the beta deriva- 
tive of the lattice classical action. Remarkably Fig. 6 shows that in the peak region the Tg 
Abelian vortex displays a signal higher than the Abelian monopole. This is better appre- 



ciated if we look at the disorder parameter Eq. ( 1.291 ) (see Fig. 7) , where In for Abelian 
vortices has a sizeable faster decrease. In the study of Abelian monopole condensation we 
found [7] that the color direction 8 is slightly preferred with respect to the direction 3 in the 
color space. It is worthwhile to see if this result holds also for Abelian vortex condensation. 
In Fig. 8 we compare the free energy derivative for the and Tg Abelian vortices for the 
lattice with = 4 and Lg = 32. As expected the Tg Abelian vortex displays, in the peak 
region, a signal about a factor two higher. 

Finally in Fig. 9 and in Fig. 10 we confront the T3, T^a, Tg and Tg Abelian vortices. 
We find that T3 and Tg agree within statistical errors in the whole range of /3 with T^a and 
T3 respectively. For the sake of completeness, in Fig. 11 we compare T3 and Tg monopoles, 
showing that they agree within statistical errors. 

We may conclude that, even for SU(3) gauge theory, our results strongly suggest that 
Abelian vortices play a role in the dynamics of confinement. 



5. Conclusions 



We investigated Abelian monopole and Abelian vortex condensation in U(l), SU(2) and 
SU(3) lattice gauge theories. 

For U(l) pure lattice gauge theory we found that, in the confined phase, the vacuum can 
be interpreted as a coherent condensate of magnetic monopoles. On the other hand, we do 
not find convincing evidence of condensation of vortices. 

For non Abelian SU(2) and SU(3) lattice gauge theories at finite temperature, by means 
of a lattice thermal partition functional, we introduced a disorder parameter for detecting 
Abelian monopole and Abelian vortex condensation in the confined phase. The disorder 
parameter is defined by means of a lattice thermal partition functional and is invariant for 
gauge transformations of the external background (monopole or vortex) field. 

Our numerical results suggest that the disorder parameter for both Abelian monopoles 
and Abelian vortices is different from zero in the confined phase and tends to zero when 
approaching the critical coupling in the thermodynamic limit. Therefore in SU(2) and 
SU(3) Abelian vortices could play a role in the dynamics of confinement. 
In particular for SU(2) gauge theory there is also a quantitative agreement between the 
measured values of free energies for monopoles and vortices. 
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On the other hand, remarkably, in SU(3) gauge theory it turns out that the Abehan vortex 
displays a signal higher than the Abelian monopole. Moreover, for the Abelian vortices 
we find that the non perturbative vacuum reacts moderately strongly to Tg vortices with 
respect to T3 vortices. This last point is in accordance with our finding in the study of 
SU(3) Abelian monopole condensation [7]. 

In conclusion it is worthwhile to observe that our results point to a difi'erent mecha- 
nism of confinement for U(l) lattice gauge theory with respect to SU(2) and SU(3) gauge 
theories. Indeed, in the U(l) Abelian case we find that the confining vacuum behaves as a 
coherent condensate of Dirac magnetic monopoles. In SU(2) and SU(3) it seems that there 
is condensation of Abelian magnetic monopoles and Abelian vortices. So that in SU(2) 
and SU(3) gauge theories one could look at the confining vacuum as a coherent Abelian 
magnetic condensate. Even more, for the SU(3) theory, it turns out that gauge field con- 
figurations leading to Abelian magnetic flux tubes, multiple of the elementary flux, seem 
to be favorite with respect to Abelian magnetic monopoles. 

We would like to remark that it is important to perform a finite size scaling analysis 
both for SU(2) and SU(3). Indeed a finite size scaling analysis will allow us to determine 
the critical behavior of the disorder parameter in the thermodynamic limit. 

We stress finally, that it should be interesting to extend our method to study center 
vortex condensation in SU(2) and SU(3) lattice gauge theories. First results in this direction 
have been reported in Ref. [52]. 
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